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V^ Abstract 

> 

^i^ Constitutive models for concrete based on the microplane concept have 

^^ repeatedly proven their abihty to well-reproduce its non-linear response 

^O on material as well as structural scales. The major obstacle to a rou- 

■^ tine application of this class of models is, however, the calibration of 

microplane-related constants from macroscopic data. The goal of this 
paper is two-fold: (i) to introduce the basic ingredients of a robust inverse 
procedure for the determination of dominant parameters of the M4 model 
proposed by Bazant and co-workers in [4J based on cascade Artificial Neu- 
ral Networks trained by Evolutionary Algorithm and (ii) to validate the 
C^ proposed methodology against a representative set of experimental data. 

The obtained results demonstrate that the soft computing-based method 
is capable of delivering the searched response with an accuracy comparable 
to the values obtained by expert users. 
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1 Introduction 



Despite long-term use of concrete in civil engineering industry, development of an 
exhaustive constitutive model for concrete still remains in the focus of engineering 
materials science. The difficulty in constitutive modeling is directly linked to 
the complexity of the material itself: the quasi-brittle response of concrete as 
seen on macroscale is a consequence of non-linear interactions of its numerous 
constituents over a wide range of lengthscales. As a result, a formulation of 
a concise model reproducing experimental data while satisfying the constraints 
of continuum thermodynamics is far from being trivial. 

One of the most successful modeling approaches stems from the microplane 
concept, in which the response of a macroscopic materials points results from 
contributions of planes of all possible orientations, see e.g. [HI Chapter 25]. 
The potential of this modeling paradigm was perhaps best demonstrated by 
the M4 model for concrete introduced by Bazant and co-workers jl], which has 
repeatedly shown its capacity to realistically reproduce responses of complex 
three-dimensional structures under general loading conditions. In addition, an 
efficient parallelization strategy [26] or model adaptivity [31] were proposed to 
compensate for an increased computational cost of the microplane models with 
respect to traditional approaches. 

The major strength of the microplane concept - the implicit format of the 
macroscopic constitutive law - is however closely linked to its major obstacle: 
the model constants are directly linked to the microplane level and as such are 
difficult to interpret and identify from experiments performed at the macroscopic 
level. In the particular case of the M4 model for concrete, Caner and Bazant 
proposed in [6] a heuristic sequence of experiments to calibrate adjustable pa- 
rameters of the model. Nevertheless, the procedure is still based on a hand- 
fitting procedure, thereby requiring an expert user or a tedious trial-and-error 
calibration procedure. To relax this constraint, hereafter we introduce a robust 
automated procedure for the calibration of the M4 model constants. 

An independent certificate to the non-trivial character of the automated M4 
model calibration is provided by difficulties associated with the selection of an 
appropriate class of optimization algorithm to solve the resulting optimization 
problem based on the minimization of the objective function. Following the 
well-established classification of the inverse algorithms introduced by Mahnken 
in [20] , the most efficient solution strategies are typically based on gradient-based 
optimization methods. Although several successful applications to quasi-brittle 
materials have been reported in the past, e.g. [2T], they cannot be applied in 
the microplane setting since a closed-form tangent is not available, resulting in 
a non-smooth objective function. 

The requirements of smoothness can be relaxed by resorting to the evolu- 
tionary algorithms. Although these methods were successfully applied to inverse 
analysis of a variety of engineering constitutive models, see e.g. [15] and references 



therein, their extremely high demand for a number of function calls combined 
with computational cost of the microplane model prohibits their use even in the 
parallel environment, see [19] . 

To reduce the high computational cost of the resulting problem, one may rec- 
oncile oneself to an approximate solution by resorting to an inexpensive approx- 
imation of the quantity of interest. At this step, network-type approximations 
are frequently employed due to their simplicity and extremely high degree of 
flexibility. In general, two complementary strategies to an approximate identifi- 
cation problem can be distinguished [15] : (i) the forward mode, where either the 
systems' response or the objective function itself is to be replaced and (ii) the in- 
verse mode where an approximate map between measurable quantities and model 
parameters is to be established. The former model was, e.g., successfully applied 
to complex geotechnical problems [29j or to an identification of parameters for 
damage-plasticity models from large-scale tests by Kucerova et. al [16]. Exten- 
sions towards the microplane case, however, suffer from the fact that an accurate 
representation of the load-displacement curve or the objective function itself still 
requires a significant number of degrees of freedom, which, when to be deter- 
mined reliably, leads to an excessive number of simulations. Moreover, a high 
accuracy of approximations is usually needed to obtain the identified parameters 
with an acceptable accuracy [16]. Therefore, high computational costs make it 
difficult to successfully implement this strategy for the microplane identification 
problem. 

The inverse identification mode, on the other hand, is substantially more 
computationally efficient as is being widely used in a wide range of material 
calibration problems, see again [15] and references therein. The critical issue in 
this case remains the inherent ill-conditioning of the problem, especially when 
determining a parameter with a low sensitivity for the executed experiment. 
This fact was exactly the reason for the failure of our early applications of the 
inverse mode to the microplane identification problem, see [S] for a particular 
illustration of this issue. Fortunately, efficient small-size procedures are currently 
available to accurately estimate the sensitivity of a particular parameter and, as 
demonstrated by our recent work [13j, present a natural choice to identify the 
model parameters. 

In the current work, we introduce a well-defined procedure for the identifi- 
cation of material parameters for M4 microplane model. The procedure itself 
is based on work by Novak and Lehky [21] and has been successfully verified 
against artificially generated sets of data in [18]. Here, we present an indepen- 
dent validation and comparison with expert user guesses. 

The remainder of the paper is organized as follows. The next section is de- 
voted to the short introduction and a description of our implementation of M4 
microplane model. The inverse mode of parameters' estimation is presented 
next with an emphasis put on artificial neural networks and their training. The 
main part is describing the methodology behind the identification of all perti- 



nent M4 model parameters. To get reliable results, two "tricks" are used. The 
first one is usage of cascade neural networks |32] where some already identified 
parameters serve as inputs to later stages of identification. The second one is 
the combination of two models for closely interacting parameters. Their com- 
bination creates nonlinear equation which is easily solvable. The main part is 
divided according to M4 authors to three sections describing three tests needed 
to estimate M4 parameters, namely the uniaxial, hydrostatic and triaxial tests, 
respectively. In the closing section, the presented text is accompanied by the 
analysis of the computational demands which present the main drawbacks of the 
presented methodology. 



2 COMPUTATIONAL MODEL 

In contrary to traditional approaches to constitutive modeling, which build on 
description via second-order strain and stress tensors at individual points in 
the (x, y, z) coordinate system, the microplane approach builds the descriptions 
on planes of arbitrary spatial orientations - so-called microplanes, related to 
a macroscopic point, see Figure [TJ This allows to formulate constitutive equa- 
tions in terms of stress and strain vectors in the coordinate system (1, m, n) 
associated with a microplane oriented by a normal vector n. The general pro- 
cedure of evaluation of a strain-driven microplane model response for a given 
"macroscopic" strain tensor £(x) can be described as follows: (i) for a given mi- 
croplane orientation n normal "macroscopic" strain tensor £(x) is projected onto 
the normal "microstrain" vector £:(n) and the shear microstrains £:(m) and e{l), 
(ii) the normal and shear microstresses cr(n),cr(m) and cr(l) are evaluated us- 
ing microplane constitutive relations, (iii) the "macroscopic" stress tensor cr(x) 
is reconstructed from the microscopic ones using the principle of virtual work, 
see, e.g., [HI Chapter 25] for more details. In the particular implementation, 
28 microplanes with a pre-defined orientation on the unit hemisphere is used to 
evaluate the response of the model. 




Figure 1: Concept of microplane modeling 
To close the microplane model description, the appropriate microplane con- 



stitutive relation must be provided to realistically describe material behavior. 
The model examined in the current work is the microplane model M4 [3]. The 
model uses volumetric-deviatoric split of the normal components of the stress 
and strain vectors, treats independently shear components of a microplane and 
introduces the concept of "boundary curves" to limit unrealistically high values 
predicted by earlier version of the model. As a result, the strain-to-stress map 
£(x) I— )■ cr(x) is no longer smooth, which complicates the formulation of consis- 
tent tangent stiffness matrix [26] and, subsequently, gradient-based approaches 
to material model parameters identification. 

In overall, the microplane model M4 needs eight parameters to describe a cer- 
tain type of concrete, namely: Young's modulus E, Poisson's ratio z/, and other 
five parameters (fci, k2, /cs, ^4, C2o]q which do not have a simple physical inter- 
pretation, and therefore it is difficult to determine their values from experiments. 
The only information available in the open literature are the bounds shown in 
the Table [H 



Parameter 




Bounds 


E 


G 


(20.0,50.0) GPa 


V 


e 


(0.1,0.3) 


h 


G 


(0.00008,0.00025) 


k2 


G 


(100.0,1000.0) 


h 


G 


(5.0,15.0) 


h 


G 


(30.0,200.0) 


C20 


G 


(0.2,5.0) 



Table 1: Bounds for the microplane model parameters 

In the present work, the computational model of a structure is provided by 
the object-oriented C++ finite element code OOFEM 1.5 [281 E?]. Spatial dis- 
cretization is performed using linear brick elements with eight integration points. 
The arc-length method with elastic stiffness matrix is used to determine the load- 
displacement curve related to the analyzed experiment. 



3 PARAMETER IDENTIFICATION METHOD- 
OLOGY 

The problem of model parameters identification starts by the design and realiza- 
tion of a suitable experiment E^x^) = y^, which connects material properties of 
our interest x^ with some observable quantities y^. Then the appropriate nu- 
merical model M(x^) = y*^ approximating the experiment in an efficient and 



^The addition of the constants Cig - C20 into M4 formulation and addition of C20 into the 
parameters' set is discussed in ^23]. The inclusion of C20 is also confirmed by the verification 



phase in Section 4.1 of this contribution. 



accurate way has to be chosen. In the next step, the model parameters x^ need 
to be cahbrated in order to minimize the difference between the model outputs 
y^-'^ and observed data y^ . 

As it was mentioned in introduction, many different identification methodolo- 
gies are available in literature. Nevertheless, only a limited number of them were 
reliably verified. As the authors of [5] introduced terms verification and valida- 
tion for numerical models of physical events, we will follow here the definition of 
those terms in the field of parameters identification presented in |15j : 



Verification: The process of determining whether the identification method 
is able to re-find the model parameters x*^ from the outputs y'^'^^ of the 
reference simulation done for any choice of original inputs x*"*^-^. 

Validation: The process of determining whether the identification method 
is able to find the model parameters x*^ corresponding to the experimental 
outputs y-^. 

In this paper, we focus on an inverse mode of an identification strategy, which 
was already verified for parameters identification on microplane model p^ and 
we would like to extend this work by validation of experimental data. Therefore, 
we just briefly review the main principles of the applied methodology and an 
interested reader can find more information in [18]. 

The methodology assumes an existence of an inverse relationship between 
model outputs y^"'^ and model inputs x*^, i.e. there is an inverse model M^^^ 
associated to the model M, which fulfils x*^ = M^^^(y*^). Generally, this 
inverse model does not need to exist. Nevertheless, we assume that the inverse 
model can be found sufficiently precise on some closed subset of the parameter 
domain. Next, we will limit our attention to an approximation of the inverse 
relationship, not its exact description. A quality of this approximation is easy 
to measure since a pair x, y obtained using inverse model M^^^(-) should also 
fulfill the forward model relation M(-). Final usage of this methodology is trivial 
because a desired value x*^ can be simply obtained as x^^ = M^^^{y^). 

The main advantage is clear. If an inverse relationship is established, then 
the retrieval of desired inputs is a matter of seconds even if executed repeatedly. 
This can be utilized for frequent identification of one model. On the contrary, 
the main disadvantages are an exhausting search for the inverse relationship, 
the existence problems for the whole search domain and inability to solve the 
problem of several global optima. 

The presented identification strategy is based on an approximation of inverse 
relation using an artificial layered neural network (ANN), which is known as 
a very general and robust approximation tool. Individual steps of the identifica- 
tion strategy involve: 

Step 1 Setup of an experimental test used for the identification procedure and 
storing the observed quantities y^. 



Step 2 Formulation of an appropriate computational model M. Input data to 
the model coincide with the parameters to be identified. 

Step 3 Randomization of input parameters. Input data are typically assumed 
to be random variables uniformly distributed on a given interval, which 
roughly lead to equally distributed precision of an ANN over the parameter 
domain. Of course, any other distribution for parameters is admissible to 
improve the ANN's accuracy e.g. in the vicinity of particular parameter 



values. A representative set of input vectors ^Train 



"■1 5 ^2 ) • • • ) ^ntr 



is carefully chosen for ANN training following the design of experiments 
methodology [22]. Another set of input vectors X;^^^ = x*^, x^, . . . , x^{g 
is randomly chosen for ANN testing, ntr and nte denote the number o 
training and testing samples, respectively. 

Step 4 Training and testing data sets preparation. The computational model 
M is applied to simulate the experiment E for all training and testing 
input vectors in order to obtain corresponding output vectors y^Train 

Yl ) y2 ) • • • ) y-ntr 



M ,.Af ,.M ^^^ ^M 



Test 



Yl ) y2 5 • • • ) Ynte 



respectively, where 



yf = M(xf ). 



Step 5 Global sensitivity analysis using the training simulations. This provides 
us with relevant model parameters which can be reliably identified from the 
computational simulation. Usually this step is performed by calculation the 
correlation between inputs X^^^j„ and outputs Yyf^^j„. 

Step 6 Definition of the topology of an ANN used for the identification proce- 
dure. 

Step 7 Training of the ANN, i.e. developing M^^^ . Usually an optimization 
algorithm is needed to appropriately setup values of synaptic weights w of 
the ANN by minimizing the error function E{w) = ||(M^^^(Y^^^j^,w) — 

vM II 
TrainW ' 

Step 8 Verification /of the ANN with respect to the computational model. This 
step is usually performed by comparing the ANN's prediction of model 
input parameters ^xist = ^^^^C^TestJ^)^ with the original one ^xist ^'^^ 
unseen testing (or reference) data. 

Step 9 Verification II of the ANN with respect to the computational model. In 
this step, a computational model should be evaluated for predicted values 
^Xest i^ order to obtain corresponding model outputs Y^^^^. Then the 
outputs Y^gg^ are compared with the original ones Y^g^^. This step is not 
necessary, but is utmost recommended. 

Step 10 Validation of the ANN with respect to the experiment. The trained 
ANN is evaluated for experimental data y^ in order to obtain correspond- 
ing input values x^ = M^^^(y^) of the computation model M. The model 



M is then evaluated for obtained inputs x"^ and results y^ = M(x^) are 
compared with original measured data y^. 

For clarity, the scheme of such identification procedure is displayed in Figure |2} 




Figure 2: Scheme of inverse analysis procedure 



3.1 NEURAL NETWORKS ARCHITECTURE 



Multi-layered neural networks became very popular as an universal approxima- 
tion tool [1]. The ANN is composed of a set of simple operational units - neurons. 
Each neuron is represented by an activation function, which defines the relation 
between an input and output value of the neuron. In the multi-layered neural 
network, neurons are ordered into layers, where outputs of neurons from one layer 
multiplied by appropriate synaptic weights serve as input values for a neuron in 
the next layer, see e.g. |S] for more details. 

When looking for a suitable neural network, two non-trivial tasks must be 
solved [33]: (a) a choice of an appropriate architecture and (b) a search for 
optimal values of synaptic weights, i.e. so-called training of a neural network. In 
our applications, we follow the premise that only three layers with enough hidden 
neurons are sufficient to describe any relation [7], therefore only three layers are 
used in our computations. Another premise used is that it is easier to train 
several networks with one output value than one complex network with several 
outputs [12j. Hence, one neural network is trained for each model parameter and 
the inverse model is then composed by a set of simple ANNs. 



Another trick used in our computations is a cascade neural network. Here, 
a sequential way of identification of the individual inputs Xi is utilized to in- 
crease the prediction abilities of the ANN-based methodology by reducing the 
complexity of the approximated relationship. The predictions of some parame- 
ters identified in the preceding steps are assumed as known and serve as inputs 
during the development of inverse models in following steps. Several applications 
of this methodology to parameters identification are presented e.g. in [32j . 

Finally, a number of input neurons that corresponds to the number of relevant 
observable quantities should be selected. Since in many engineering problems, 
measurements are represented by single or several curves usually discretized into 
discrete points, a number of observed quantities is often very large. Therefore, 
a global sensitivity analysis is applied to choose only quantities important for 
a particular model parameter to be identified |2l]. A number of neurons in 
hidden layer is determined by consecutively increasing their number taking the 
over-training and under-training issues into account. Particularly, a goal is to find 
a trade-off between errors in prediction for training and testing data, respectively. 



3.2 NEURAL NETWORK TRAINING 



Once the architecture of the ANN is chosen, the training process can start. The 
most famous training algorithm still remains a relatively simple backpropagation, 
especially because of its simple interpretation and implementation. As the name 
suggests, the error in prediction of the output layer is propagated back to previous 
layers and according to an error's value, synaptic weights are updated. From the 
optimization point of view, it is a gradient-based algorithm, which suffers from 
premature convergence. 

Therefore, GRADE evolutionary algorithm extended by a niching strategy 
CERAF is used here as a more robust optimization algorithm. This algorithm 
was derived from its predecessor SADE algorithm, which provided better results 
in ANN training against the back-propagation as was presented in [S]. GRADE 
is a population-based real-coded algorithm consisting of three genetic operators, 
which are applied to a set of solutions. Particularly, it involves mutation, dif- 
ferential cross-over and inverse tournament selection. CERAF is a multistart 
strategy enhanced with memory that increases the ability of a population-based 
algorithm to escape from local extremes. An interested reader can find more 
details about both methods in [15] or [171 [Tl], where C++ and MATLAB imple- 
mentations are available. For information on other possible training algorithms 
see e.g. [12] or a review [33] . 



4 EXPERIMENTAL VALIDATION 

Following the heuristic calibration procedure suggested in |6], we examine three 
specific experimental tests: (i) uniaxial compression, (ii) hydrostatic test and 
(iii) triaxial test. Advantage of these tests is their simplicity and availability in 
most experimental facilities. Moreover, authors in [6] claim that these experi- 
ments are sufficient to determine all parameters of the microplane model M4. 
The results presented in this section can be understood as a verification of this 
claim. 

In our previous contribution [18], we have shown that the proposed methodol- 
ogy is able to identify all parameters from computer-simulated curves. To demon- 
strate the applicability of the proposed procedure, a real simulation should be 
examined. Since the measurements from different loading tests are obtained for 
different concretes, this section does not represent a validation of proposed iden- 
tification strategy in general, but only the validation of application of particular 
inverse models. 

Each neural network is obtained by the procedure described previously. More 
precisely, a fully connected three-layer perceptron with bias neurons is applied to 
map discrete values from stress-strain diagrams to microplane model parameters. 
Log-sigmoid functions are considered as activation functions in all neurons. In all 
cases studied in this work, 60 training samples are generated by Latin Hypercube 
Sampling method and optimized by Simulated Annealing in order to minimize 
the correlation among all samples p^. Those two methods are implemented in 
FREET software [25J that has been used. Next, 10 testing samples are obtained 
for random parameters from given bounds. Relevant values selected from stress- 
strain diagram as ANN's inputs are chosen by hand with respect to the results 
of global sensitivity analysis. Here, the Pearson product-moment correlation 
coefficient [10] is evaluated between the discrete values of stresses (or strains) 
and corresponding values of microplane model parameters. For the ANN training 
the GRADE algorithm supported by CERAF strategy is used and calculation is 
stopped after 1,000,000 iterations of the algorithm. Finally, all three following 
compression tests are performed on the cylinders with a radius equal to 75 mm 
and the height of 300 mm. 



4.1 UNIAXIAL COMPRESSION TEST 

Experimental data are taken from [6], Figure la. Note that the original source 
comes from [30]. A suite of 70 training and testing simulations and the resulting 
bundle of curves together with the experimental data is shown in Figure |3^. Only 
for illustrative purposes, the evolution of global sensitivity during the loading 
process (horizontal axis) is depicted in Figure 111 and is omitted for other two tests. 
The results indicate that the most sensitive parameters are Young's modulus E, 
the coefficient fci, Poisson's ratio ly (especially for the initial stages of loading) 
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0.004 0.006 

e[-] 



(a) 



(b) 



Figure 3: Bundle of curves of uniaxial compression test simulation with experi- 
mental data for (a) all parameters and (b) only constant C20 identification. 

and, for the later stages of loading, the coefficient C2o- Therefore, one can expect 
that only these parameters can be reliably identified from this test. 
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Figure 4: Sensitivity evolution for uniaxial compression test. 



Moreover, the impact of individual parameters on a position of a peak of stress- 
strain curves is computed. The results of a sensitivity analysis using Pearson's 
product moment correlation coefficient of peak coordinates [epeak,Cpeak] are plot- 
ted in Fig |5| Results indicate particularly strong influence of the ki parameter, 
which promise its reliable determination. 

Based on the results of sensitivity analysis, the neural network training can 
be performed using a nested strategy of cascade neural networks. First, Young's 
modulus E with sensitivity f^ 1 in the initial stage is easily identifled. To this 
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Figure 5: Pearson's coefficient as sensitivity measure of individual parameters to 
peak coordinates [epeaki^^peak] of stress-strain curves 

end, the following ANN's inputs are chosen: the values of stresses cr[e], where e is 
a particular level of deformation with extremal Pearson's correlation coefficient, 
here the starting of the loading curve, particularly e = 0.0005; 0.001 and 0.0015, 
respectively. The hidden layer contains two neurons only; the output layer con- 
sists of one neuron corresponding to the predicted value of Young's modulus E, 
see also notation used in Table [2j Similar approach is used for ki parameter with 
the architecture described in Table [2] Note the usage of Young's modulus E 
estimation from the previous ANN. Errors for both ANN's predictions are then 
listed in Table [3l 



Parameter Topology 



Inputs 



E 3 + 2 + 1 o"[o.ooo5]) 0" [0.001]) o" [0.0015] 

kl 5 + 3 + 1 (T[0.0025] , q'[0.009] , (-peak, C^peak, prcdictioU of E 

Table 2: Topology of neural networks for uniaxial compression test. 

Next, neural networks were applied on measured data and following values 
were obtained: 



E 

h 



36057MPa 
0.000196 



Since the lateral deformation for measured data are missing. Poison's ratio 
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Parameter 



Training data Testing data 

Maximal error Average error Maximal error Average error 



E 



1.97 1.01 1.50 0.81 

2.68 1.57 2.63 1.50 



Table 3: Errors in ANN's prediction relative to the definition interval of the 
parameters in [%]. 

cannot be identified. Its value is chosen according to microplane authors' rec- 
ommendations V = 0.2. Concerning other parameters, according to sensitivity 
analysis, other parameters should not be important. Therefore their values are 
set to mean values of predefined intervals: 



k2 -- 


= 550 


fcs = 


= 10 


k, -- 


= 115 



Finally, recommended value of a constant C20 by M4 authors is C20 = 1-0. 
In Figure |6} the measured data are compared to the simulation obtained for 
predicted parameters and to the simulation obtained by M4 authors using hand- 
fitting method. 

Following our experience [23j, it is interesting to predict also the value of 
constant C2o- Therefore, new 70 simulations (60 for training and 10 for testing 
purposes, respectively) were performed with fixed values of E, v and k\ param- 
eters. The resulting bundle of curves can be compared to measured data, see 
Figure |3|d. The same discretization was respected. The topology and inputs for 
the neural network are presented in Table |4] and errors in ANN's predictions are 
listed in Table [H 

Parameter Topology Inputs 

C20 4 + 2 + 1 O" [0.003] ; O" [0.004] ; O" [0.006] ; ^[0.008] 

Table 4: Topology of neural networks for constant C20 identification. 



Parameter 

C20 



Training data Testing data 

Maximal error Average error Maximal error Average error 
11.18 6i)2 11.27 (U5 



Table 5: Errors in ANN's prediction of constant C20 relative to the definition 
interval of the parameters in [%] . 

Neural networks were applied on measured data and following value was ob- 
tained: C20 = 0.72785. 

The final comparison of measured data and predicted simulations is shown 
in Figure |6j To show some more objective comparison of presented simulations. 
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Figure 6: Comparison of measured data and predicted simulations for uniaxial 
compression test. 

the error between measured data and simulated curves can be summarized over 
discrete points corresponding to measured data, i.e. 



E 



(1) 



where a = [. . . , a^, . . .] corresponds to stresses measured in discrete points i and 
(T = [. . . , cTj, . . .] are corresponding simulated values. Values of the error defined 
by Equation [T] are shown in Table [6] 



Prediction 



Error E 



expert estimation jSj 8.30 

ANN + C20 = 1 12.88 

ANN + C20 = 0.72785 6.97 

Table 6: Comparison of errors of predicted simulations. 



4.2 HYDROSTATIC COMPRESSION TEST 

The next independent test used for the identification problem is the hydrostatic 
compression test, where a concrete cylinder is subjected to an increasing uniform 
pressure. Experimental data come again from [6J. These data were obtained by 
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authors Green and Swanson [9]. The stress-strain diagram represents the rela- 
tion of hydrostatic pressure a and axial deformation e. The results from simula- 
tion performed for microplane model parameters obtained by expert estimation 
are there also available in comparison with measured data. Since we suppose 
that values of Young's modulus, Poisson's ratio and ki parameter can be reli- 
ably obtained from identification of uniaxial compression test (these results are 
nevertheless not available for the concrete observed here), the values of these 
parameters are taken directly from the article [6J. Then, the goal is to identify 
values of parameters k^ and k^\ parameter ^2 attains zero sensitivity here. The 
bundle of resulting stress-strain diagrams for training and testing sets can be 
compared with measured data in Figure [7j 



500 



^ 300 



— simulations 
'•-^ measured data 




0.05 0.1 0.15 0.2 0.25 



Figure 7: Comparison of measured data and results of 70 simulations of hydro- 
static compression test. 



Parameter 



k. 



ANN'S layout 



3-2-1 



Input values 



^3; ^peaki e[85.5],-u 



Table 7: Neural network architecture for ^4 parameter identification; e[85.5],« is the 
deformation corresponding to hydrostatic stress a = 85.5 MPa in the unloading 
phase. 

During the verification process, we have not been able to properly identified k^ 
and ^4 parameters. Since all other parameters have very small or zero correlation, 
we suppose that the nonlinearity in the parameter k^ identification is caused by 
k^ parameter and vice-versa. In other words, these difficulties can be caused 
by some level of correlation between the parameters k^ and ^4. Therefore, to 
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Table 8: Des 


Topology Inputs 


ANNi 4 + 2 + 1 /C4, epeafe, ^yield, €[214],/ 
ANN2 5 + 2 + 1 A;4, epeak, (-yield, €[137],;, €[308],/ 


cription of two neural networks trained to predict A;3 parameter 


Parameter 


Training data 

Average error Maximal error 


Testing data 

Average error Maximal error 




1.40 2.59 
1.51 2.52 


1.71 3.07 
1.21 2.13 



Table 9: Error in ANN's predictions relative to the definition interval of the 
parameters in [%]. 

eliminate unknown correlation between parameters k^ and ^4, their values are 
used also as inputs into ANN's. For k^ parameter, an ANN was then constructed, 
see its architecture in Table [7l 

The inputs to the ANN for prediction of k^ parameter were chosen in two 
different ways in order to show an impact of selection particular points of the 
post-elastic loading part of the stress-strain diagram. Two neural networks were 
trained with the deformation value (i) in the middle of loading phase with stress 
a = 214 MPa and (ii) in its thirds: a = 137 MPa and a = 308 MPa leading 
to topologies described in Table [8j Hence, the second ANN puts more impact 
on the post-elastic part of the diagram than the first one. Both neural networks 
have among the inputs the value of k^ parameter and values of deformation at 
the end of loading and at the end of the elastic part of the diagram. 

Errors in predictions for training and testing data are noted in Table [9j For 
the sake of brevity only errors corresponding to the first ANN trained for ^3 
parameter are presented, the errors of the second network were negligibly higher. 

Two neural networks obtained to predict k^ and k^ parameters with inputs 
taken from measured data represent system of two non-linear equations, which 
can be solved graphically. Relation given by a neural network for /C4 prediction 
with all inputs fixed to values obtained from measured diagram except the input 
value of ^3 parameter, which states a variable, is shown in Figure |8j The relation 
is defined in similar manner by the first neural network trained to predict k^ 
parameter. The intersection of presented curves defines the predicted values of 
ks and k^ parameters corresponding to measured data, in this case ^3 = 12.34 
and ^4 = 98.19. Same Figure also shows an equivalent relation for k^ parameter 
obtained by the second neural network. Coordinates of the intersection are k^ = 
11.20 and k^ = 91.85 in this case. 

In both cases, the predicted values differs from values stated by authors in [H], 
where ^3 = 9 and k^ = 82. The comparison of measured data with simulated 
diagrams obtained for parameters given in [6] and for two couples of parameters 
predicted by neural networks is depicted in Figure [9] 
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k^=ANN(k3) 
k3 = ANN^(kJ 
k3 = ANN2(kJ 



40 60 80 100 120 140 160 180 200 



Figure 8: Relations of k^ and ^4 parameters for measured data: red continuous 
curve represents ANN trained to predict ^4 parameter, black dash-and-dot curve 
corresponds to ANN trained to predict k^ parameter with four inputs (ANNi) 
and blue dash curve to five inputs (ANN2). 
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Figure 9: Comparison of measured data and simulated diagrams of hydrostatic 
compression test for predicted parameters. 



It is not so easy to judge which predicted simulation really better corresponds 
to measured data. A simulation presented in [6j relatively well fits the measured 
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data at the end of elasticity and at the end of loading. Nevertheless, there is 
a significant error in the middle of the loading part of the diagram. A simula- 
tion corresponding to the first prediction of k^ parameter relatively well fits the 
measured data in the middle of the loading diagram and is worse near the end of 
elasticity and the end of the loading. A simulation corresponding to the second 
prediction of k^ parameter can represent a compromise between the two previous 
simulations. To show more objective comparison of presented simulations, the 
error between measured data and simulated curves can be evaluated similarly 
to Eq. ([I| just in terms of deformation instead of stresses. The values of the 



resulting error are presented in Table 10 



Prediction Error E 

expert estimation [6] 0.0257 
ANNl 0.0163 

ANN2 0.0161 

Table 10: Comparison of errors of predicted simulations. 



Prom the comparison presented in Table pi)] it is clearly visible that simulations 
performed for parameters predicted by both neural networks fit the measured 
data better than the simulation done for parameters obtained by expert estima- 
tion presented in |6]. Moreover, it is also visible that simulations predicted by 
neural networks are somehow handicapped. Implementation of the microplane 
model M4 in the OOFEM software does not properly describes non-linear be- 
havior in the unloading part of the diagram that leads to a higher error. 

4.3 TRIAXIAL COMPRESSION TEST 

The last experiment, used for the purpose of parameter identification, is a triaxial 
compression test. To this end, a specimen is subjected to the hydrostatic pres- 
sure an- After the peak value of an is reached, the axial stress is proportionally 
increased. The "excess" axial strain e = St — Sh, where St and Sh denote the to- 
tal and hydrostatic axial strain, is measured as a function of the overall stress a. 
Similarly to the hydrostatic compression test, also for triaxial compression test 
we use measured data from [B] . These data were obtained by Balmer ^ . Again, 
the presented measured data are accompanied by the simulation performed by 
the authors of [6] for parameters values established by the expert estimation. 
The triaxial compression test is supposed to be the last experiment needed to 
identify parameters, which cannot be identified by uniaxial or hydrostatic com- 
pression tests, i.e. ^2 parameter. Therefore, Young's modulus, Poisson's ratio, 
/ci, /c3 and k^ parameters are supposed to be known and their values are taken 
directly from [6|. 

Pive different measurements for the triaxial compression test are available 
in [6] corresponding to five different levels of a hydrostatic pressure an applied 
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to specimens (an e {34.5,68.9, 103.4, 137.9, 172.4} MPa). The bundle of 70 re- 



sulting stress-strain diagrams can be compared with measured data in Figure 10 
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Figure 10: Comparison of measured data and results of 70 simulation of triaxial 
compression test. 



In Figure 10 it is clearly visible, that measured data are remote from all 
simulated curves. This can be caused by wrong limits for k^ parameter, i.e. 
k2 G (100; 1000), however an increase of the upper bound by a factor of two did 
not bring substantial shift towards experimental data [15]. The discrepancy can 
be also caused by other parameters, i.e. Young's modulus. Poison's ratio, fci, k-^ 
or k/^ parameter. Nevertheless, the goal of this Section is not to identify these 
parameters from triaxial compression test and therefore, we will test the ability 
of a neural network to cope with this problem. 



Parameter 



u 



ANN'S layout 



Input values 



3-2-1 



Cpeafc, tX [0.0128]) '^[0.0308] 



Table 11: Neural network architecture for k2 parameter identification. 



A neural network with an architecture presented in Table 11 was trained on 
simulated data and then applied to predict a value of /c2 parameter for measured 
data. Note that only data from the lowest loading level were used. Errors of 



Ann's predictions on training and testing samples are written in Table [12 

The prediction of the neural network for measured data is ^2 = 1193. It is 
not surprising that the neural network needs to extrapolate and its prediction 
exceeds the limit given for ^2 parameter. Although layered neural networks are 
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Parameter 



Training data Testing data 

Average error Maximal error Average error Maximal error 



k. 



2.63 



5.82 



2.37 



6.23 



Table 12: Error in ANN's predictions relative to the definition interval of the k2 
parameter in [%]. 



known to be good in approximations, they are week in extrapolation. Moreover, 
there is a question, whether it is possible to find appropriate value of ^2 param- 
eter to fit measured data, since it is probable that a more important error is 



hidden elsewhere. Finally, Figure 11 shows the comparison of measured data, 
a simulation given in [B] and a simulation for the parameter value predicted by 
the neural network. 
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X measured data 
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Figure 11: Comparison of measured data and simulated diagrams of triaxial 
compression test for predicted parameters. 



Similarly to previous sections, errors between measured data and simulated 
curves can be calculated. Resulting values for all five experiments corresponding 



to different levels of hydrostatic pressure an are shown in Table 13 



Errors in Table 13 are very similar for both simulations. It is possible to 
conclude that even for remote input data the neural network is able to predict 
a reasonable value for ^2 parameter and resulting simulation is comparable with 
that one obtained by the expert estimation. 



20 



an [MPa 


expert 


estimation 


6 


ANN'S 


predictions 


34.5 




64 






43 


68.9 




87 






60 


103.4 




73 






45 


137.9 




70 






39 


172.4 




43 






60 



Table 13: Comparison of errors of predicted simulations for triaxial compression 

test. 

5 Conclusions 

In this contribution, an example of the engineering problem, which is difficult to 
be solved by traditional procedures, was solved using soft computing methods. 
Particularly, cascade neural networks were used to estimate required microplane 
material model parameters in a sequential way. As the training procedure, the 
evolutionary-based method GRADE extended by CERAF strategy was used. 
A number of needed simulations is reduced by the application of the Latin Hy- 
percube Sampling method accompanied by the optimization by Simulated An- 
nealing. The global sensitivity analysis shows not only the influence of individual 
parameters but also approximately predicts the errors produced by neural net- 
works. 



Parameter 


Test 


ANN'S topology 


ANN'S inputs 


E 


Uniaxial compression 


3 + 2 + 1 


(^z,l, C^z,2, 0"z,3 




V 


Uniaxial compression 


4 + 3 + 1 


O'xi, 0-X2, E, ki 




ki 


Uniaxial compression 


4 + 2 + 1 


<^x,2j 0'z,peak, ^z,peak, E 




fe 


Triaxial loading 


3 + 2 + 1 


O'peak, Cr29, 0"lOO 




ks 


Hydrostatic loading 


5 + 2 + 1 


"^4; ^yield) ^load,2) ^load,bi 


^peak 


ki 


Hydrostatic loading 


3 + 2 + 1 


"'S) ^peak) ^unload,A 




C20 


Uniaxial compression 


3 + 2 + 1 


O'peak, <761, 0-81 





Table 14: Final status of M4 identification project 



Results, see Table 14, confirm the claims made by authors [6] of the mi- 



croplane M4 model on individual parameters fitting. However, the validation of 
the complete process cannot be done, since the experimental data for all three 
loading tests performed on one concrete are not available and probably still do 
not exist. Therefore, the validation demonstrated in this paper is done only for 
particular identification steps. 

The rather severe disadvantage of the microplane model, and also of the pro- 
posed methodology, is an extreme demand of computational timq^ However, 
especially the uniaxial test dramatically depends on how far the computation 



^Note that since the task is symmetric, only 1/8 of the cyhnder is actually needed for FEM 
analysis which shortens the needed time. 
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Figure 12: Time required to access particular point on the uniaxial load- 
deflection curve. 



will be processed. If we inspect Figure 12, every 0.001 move on the softening 
branch costs in average one hour on a single processor PC with the Pentium IV 
3400 MHz processor and 3 GB RAM. However there is a problem that individual 
compositions lead to big differences in times. This is illustrated in Figure 13 
where the scatter of times is shown in a form of histograms. The averages are 
8 minutes, 19 minutes and one hour for hydrostatic, triaxial and uniaxial com- 
pression tests, respectively. However note, that in cases that the material is very 
soft and the load-deflection curve is long, some computations longer than 6 hours 
can appear. 

Because the identiflcation procedure consists of developing cascade neural net- 
works, all but one level inverse models should be recalculated for any new mea- 
sured data. Fortunately, the most time consuming simulations of the uniaxial 
compression test are necessary for training of the flrst three neural networks pre- 
dicting Young's modulus, Poisson's ratio and ki and can be used repeatedly for 
any new measurement. Then the proposed methodology still needs to compute 
30 uniaxial tests to properly identify C20 parameter and a set of 30 hydrostatic 
and triaxial tests to fit k^, k^ and ^2- 
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is measured at e = 0.004 to fit the figure. 
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